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In this research, we propose a method to calculate the phase delay angle of 
the reflected wave compared to the wave emitted from the interdigital 
transducer (IDT) of a passive wireless surface acoustic wave (SAW) sensor 
by combining simulation of the finite element method (FEM) and 
MATLAB. Based on the result, we consider the influence of the reflector 
position on the phase delay of the SAW sensor consisting of an IDT and 
three equally spaced reflectors. The distances between the IDT and the first 
reflector are 870 micrometers, 939.6 micrometers, and 1009.2 micrometers, 
and the distances between the reflector positions are 452.4 micrometers, 522 
micrometers, and 591.6 micrometers. The outputs of the FEM simulation, 
including delay time, amplitude loss, and the length of the reflected wave are 
put into the model in MATLAB to determine the phase delay of the response 
waves. The obtained phase delay results show that as the distance between 
the IDT and the reflectors increases, the phase delay angle increases to 4.49 


degrees, 5.39 degrees, and 5.78 degrees, respectively. 
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1. INTRODUCTION 

Surface acoustic wave (SAW) sensors are widely used in measurement and identification systems, 
and applicated in different fields such as measuring high temperatures in harsh environmental conditions 
[1]-[3], pressure [4], [5], torque [6], [7], magnetic and humidity [8] or in identification tags (ID-tag) by their 
superiority as no battery or do not require any power supply to operate, wireless, small, made on highly stable 
materials. Using passive wireless SAW sensors requires a reader to transmit radio frequency (RF) signals and 
receive reflected signals that carry measurement information. Several publications have mentioned the reader 
of the passive wireless SAW sensor [9]-[14]. The previous authors [9] built a sensor network model with four 
wireless SAWs for industrial applications, used a network analyzer to read the reflected pulses and calculated 
the result on MATLAB software, but they did not mention the algorithms on the reader. Bruckner and 
Bardong [10], have investigated the reader architecture, methods of time-domain sampling (TDS), and 
frequency-domain sampling (FDS), in which the author also analyzed and proposed the theory of hybrid 
structure TDS/FDS for the reader, but the article has not mentioned the result of the simulation or experimental 
on the device. In the fabrication process of a SAW temperature sensor presented in [11], the authors used the 
results from the finite element method (FEM) and the boundary element method (BEM) simulation to fabricate 
sensors, but the computational solution to read the feedback pulse has not been mentioned. 
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Some experiments on a reader have been presented. According to Krishnamurthy et al. [12], the 
transmission and reception of the reader were separated by on/off the switching at 50 Khz. The transmitted and 
received signals are simulated with a relatively long reading period, making them unsuitable for highly dynamic 
measurements. The results published in [15] have indicated a method to calculate the phase delay between the 
reflected waves and interrogation signals of the SAW delay line sensor, thereby calculating the tan @ (¢ọ is the 
phase shift angle of reflected pulses) using a model on MATLAB software. However, the data input to the 
model was random values, which had been chosen to test the model and the tang is not convenient for scale. In 
addition, the article did not consider the influence of reflectors. Research by Kozlovski et al. [16], built a reader 
to read the temperature of four SAW sensors at the same time, but the number of sensors was small and the 
manufacturing process was complicated, requiring a lot of experience to test the different properties of the 
sensor. 

For passive wireless SAW sensors, the phase shift between interrogation signals and response pulses 
contains measurement informations [17]. The reader must capture the peak of the reflected wave in order to 
calculate the phase shift corresponding to the measurement. Studies [15], [17] have shown that the phase shift 
is affected by the measurand and shown the time delay between sending the probe pulse and receiving the 
reflected pulse. In addition, the position of the reflectors also affects the delay time of the reflected wave 
[18], so it is necessary to evaluate the position of the reflector with respect to the phase shift angle. In this 
study, we look at a specific way to figure out the phase angle directly from the delay time of the reflected 
pulse, and we also look at how the position of the reflector affects the phase shift. 

Because of advantages such as easy device structure modeling, adaptability, and time-dependent 
simulation. FEM and MATLAB have been used by many researchers for design and simulation in various 
sensor fields. Kozlovski et al. [16], used MATLAB to program signal processing and temperature extraction 
in the 915 MHz correlator transceiver and sensor system. While using FEM to analyze and optimize a 
candidate ultra-wideband antenna for Ku-band wireless is mentioned in [19]. In this study, we use a 
combination of FEM simulation and MATLAB to simulate single-port SAW reflective delay line (RDL) 
structure on FEM with YZ-LINbO; piezoelectric substrate. The outputs of the FEM simulation were 
imported into blocks in MATLAB to calculate the phase difference. 


2. SIMULATION PROCESS 

SAW device configurations include SAW delay line, one-port SAW resonator, and two-port SAW 
resonator. In this paper, the authors used the SAW RDL because they are commonly used for passive 
wireless sensing systems [20]. 


2.1. Principle of passive wireless SAW 

Figure 1(a) shows the principle of the one-port SAW RDL and the phase difference of the reflected 
signals are shown in Figure 1(b). The system has two devices, the reader and the sensor. Where the sensor is 
a piezoelectric substrate on the surface put an interdigital transducer (IDT) and reflectors, and there is an 
antenna attached to the IDT, the reader also attached an antenna to transmit/receive electromagnetic signals. 
The reader emits interrogation pulses propagation to the antenna of the IDT. Here, the electrical signals are 
converted into acoustic waves and propagate on the surface of the piezoelectric substrate. When the acoustic 
waves propagate to the reflectors, a part of these waves will reflect in the IDT. The antenna will receive the 
reflected pulses and sends them back to the reader, where they should compute to determine the sensor's 
measurement values. 
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Figure 1. The one-port SAW RDL (a) the top view of the structure and (b) the phase difference between the 
response and the interrogation pulse 
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The propagation of SAWs depends on the geometry of the substrate and environmental conditions. 

In Table 1, the obtainable linear coefficients for physical effects on SAW substrates are summarized [21]. 

The response pulses from the reflectors provide the measurement informations. In Figure 1(b), œ and 9, are 

the phase shift in the first and the second response, respectively, and the interrogation pulse. The propagation 

time of acoustic wave on the substrate is: 
Li 


Ciara 


Vsaw 


(1) 


where L; is the distance from IDT to the in reflector and Vsaw is the velocity of SAWs. 
In SAW sensor applications, L; and Vsaw are changed due to variations in temperature, mechanical 
stress, and strain on the surface layer. 


Table 1. Linear coefficients for physical effects on SAW substrate materials [21] 


Physical quantity Linear coefficient 
Temperature Up to 100 ppm/K 
Pressure, stress 2 ppm/kPa 

Force 10 ppm/kN 

Mass loading 30 ppm/ug.cm? 
Voltage 1 ppm/V 

Electric field 30 ppm/V.um 


Reflectors are designed to a part of the incoming waves from the IDT will be reflected while the 
remainder will pass through and reach the next. The relative sensitivity of T for the effect parameter y [21] is: 


Sp=z. (2) 
In a linear approximation (1* order Taylor expansion) [21], we obtain: 
Ty =To(1+ S?.y) (3) 


where Tp is the delay time of the response signal from the reflector, and Ty is the delay time when measuring 
the y parameter. 
The variation AT of a time delay Tọ is: 


APTS TS, (4) 


Measuring the time delay between transmitted and received signals can be done in a variety of methods. 
Figure 1(b), the functional relationship of the temperature T and the phase shift between the peaks of the first 
and second reflections are depicted by: 


Q2-1(T) = P2-1(Tref ) [1 +TCD*(T-— Tref)] (5) 


where Tef is the reference temperature (or room temperature) and TCD is the temperature hysteresis 
coefficient of the substrate material. This phase shift is affected by environmental factors (temperature, 
pressure, humidity [22]), and the position of the reflectors [5]. Gruber et al. [15] modeled SAW on 
MATLAB to calculate tang but did not consider the influence of reflectors. 

This paper evaluates the phase shift between the feedback pulses in which we have used the results 
obtained from the FEM simulation as amplitude, wavelength, and time delay of the reflected waves to build 
the model and algorithm on MATLAB to calculate the phase shift angle. The details are presented in sections 
2.2 and 2.3. 


2.2. FEM simulation 

Figure 2(a) shows the SAW structure used in the FEM simulation. The SAW has three aluminum 
reflectors, a substrate made by YZ-LiNbOs, and an IDT. The reflector combines four electrodes. The IDT is 
bidirectional with 50 pairs of fingers, and the wavelength (A) is 39.94 um. In the applications of wireless 
passive SAW sensors, the reflectors are chosen so that the response signals are as large as possible. As the 
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thickness of the bars in the reflector increases, the amplitude of the reflected wave also increases [18]. 
However, increasing the thickness of the fingers of IDT will reduce the resonant frequency [23]. In this 
paper, we have chosen the thickness of the bars in the first reflector to be equal to the thickness of the IDT 
fingers. The thickness of the IDT finger and reflector 1 are 2 um so the ratio gives h)/A = 0.05 [24] and the 
thickness of the reflectors 2,3 is 4 um, respectively h2/A=0.1 to increase the amplitude of the 2"! and 3" 
reflected waves. 

Figure 2(b) shows the cross-sectional view of two IDT fingers, where d is the period of electrodes, 
a is the distance between the electrodes, h, is the thickness of the IDT fingers, and p is the periodicity of a 
periodic structure. Table 2 indicates the detailed parameters used in the simulation. The center frequency 
is [25]: 


f=v/à (6) 


where v is the velocity of the SAW and 4 is the wavelength. The propagation velocity of the surface wave 
depends on the shear angle of the substrate material or elastic matrix [26] as (7): 


v= = (m/s) (7) 
where 
a (8) 
IDT 
| Reflector1 Reflector? Reflector3 l p > d < 
ay az a3 a3 / b2 — 
Sc se sels v ; 


L 


(a) (b) 


Figure 2. The passive wireless SAW structure (a) the 2D structure was used in the FEM simulation and 
(b) the cross-sectional view of two IDT fingers 


With the YZ-LiNbO3 substrate, coefficients of the elastic matrix, cj;, and density are described in 
Table 3 [24]. From (7) and (8), we calculated the velocity of wave propagation on the substrate as 3994 m/s 
and wavelength A=39.94 um with the center frequency of 100 Mhz. We test the surface wave propagation 
velocity during simulation by applying a Dirac delta function to the IDT and then reading the propagating 
pulse to the first reflector to find the correct center frequency. 


Table 2. Geometric parameters used in the FEM simulation 


Parameters Structure 1 Structure 2 Structure 3 

The length of the substrate, Le (um) 5394 5463.6 5533.2 
The number finger of IDT, N 50 50 50 
The wavelength, à (um) 39.94 39.94 39.94 
The thickness of the substrate, h (um) 400 400 400 
The IDT-reflector 1 gap, a2 (um) 870 939.6 1009.2 
The reflector-reflector gap, a3 (um) 452.4 522 591.6 
The IDT-left boundary gap, a; (um) 1740 1740 1740 
The reflector 3- right boundary gap, bz (um) 1740 1740 1740 
The thickness of the IDT fingers and reflector 1, hı (um) 2, 2 2 
The thickness of the IDT fingers and reflector 2,3, h (um) 3.99 3.99 3.99 


In the structure of SAW one-port with three reflectors simulated FEM as shown in Figure 2(a), to 
reduce the simulation time in the meshing process, the substrate below the IDTs and reflectors has a thick 
mesh with 40 elements per wavelength and other regions have a sparse mesh with 10 to 20 elements per 
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wavelength. With the above meshing, the total number of elements obtained is 41516 and the solving time is 
2 us. To calculate the propagation velocity of the acoustic wave, firstly, we put the Dirac pulse in the IDT 
and read the transmitting pulse to the first reflector. After that, we take a fast Fourier transform FFT to get a 
frequency response at f=102.4 Mh- as shown in Figure 3. With this frequency, from (1.4), we can calculate 
the velocity of the transmitting wave as 3702.72 m/s. 

With the received frequency, the AC voltage applied to the IDT is 
2« cos(2m * 102.4 « 10° * t)rect(t — 50 * 10°), where t is the time and rect(x) = 1, when 0< x < 1. 
Otherwise, rect(x) = 0; T is the width of the Dirac pulse. In Figure 2, the right boundary gap and the 
reflector—the left boundary gap and the IDT are calculated to evade the reflected waves from boundaries 
since we only care about the reflected waves from these reflectors. Solution time is 2 us ensures that the 
response is received from the reflectors. 


Table 3. Elastic matrix coefficient (GPa) and density p(10? kg/m?) of YZ-LiNbO; 
Cu Ci2 C13 C33 C44 C4 p 


203 53 75 245 60 9 47 


-70 X: 1.024e+08 
Y: -73.76 


` V 


Insertion Loss (dB) 


1 1.01 1.02 1.03 1.04 1.05 
Frequency (Hz) < 108 


Figure 3. Frequency response when applying the Dirac pulse to the IDT 


2.3. The modeling on MATLAB 

The one-port SAW delay line sensor is not able to model with only one port, it has been described 
as a two-port system on MATLAB, consisting of input and output, but with the same behavior, as Figure 4. 
When putting a voltage to the IDT, due to the piezoelectric effect described in section 2.1, a SAWs will 
propagate to the reflectors 1-3. Because the reflectors are located at different positions on the piezoelectric 
substrate, the reflected waves will be obtained at different times T1, T2, and T3. Figure 4 shows the results 
from response pulses in the FEM simulation as Gain;, l;, and T; (l; are the length of the reflected waves, and 
gain; represents the amplitude loss of the reflected waves, respectively, (i = 1 + 3). These results are put in 
the block diagram of the reception and transmission model at IDT with three reflectors of SAW on 
MATLAB/Simulink. Evaluating the signal delay allows for the calculation desired measurement quantity. 
Figure 5 shows the MATLAB/Simulink model for calculating phase parameters. 

The voltage applied to the IDT is: 


u(t) = A * cos(w * t) (V) (9) 
where A is the input voltage amplitude, A = 2(V); œw is the frequency of the input signal, w = 2 * m * 
102.4 « 10°(rad/s). 

The sum of response signals from three reflectors at the output (Figure 4) is: 


dizi A * A; *cos (w * (t= T;) + 8) * JI- 7) (10) 


where A; is the amplitude loss for the reflectors after passing the it” reflector, respectively; n is the number 
of reflectors, n= 3; T; is the delay time corresponding to the it” reflector. 
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Figure 4. The results in the FEM simulation put in the block diagram of SAW on MATLAB/Simulink 
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Figure 5. The Simulink model of calculating phase parameter 


3. RESULTS AND DISCUSSION 

Figure 6 indicates the amplitude of the input wave on the IDT and reflected waves from the first, the 
second, and the third reflector in FEM simulation. These reflected waves reach the maximum amplitude at 
the times T,, T2, and T3, and have lengths l,, l}, and l}, respectively. In the simulated structure, the maximum 
amplitude of the response pulses from the first, the second and the third reflector are 0.05604 (V), 0.02687 
(V) and 0.01811 (V) at the times t_1=0.487 (us), t _2=0.806 (us), and t_3=1.132 (us), respectively. The 
times when the IDT receives the reflected waves according to the theoretical calculation are t_ 11, t 21, and 
t 31. Where t_11=(L_e/2+2*a_2)/v=0.5157(us); 1t_21=(L_e/2+(a_2+a_3+7*d)*2)/v=0.798 (us), and 
t 31=(L_e/2+(a_2+2*a_3+14*d)*2)/v=1.054(s), respectively, where L_e is the length of IDT, L_e=(4*N- 
1)*d (um). The times to reach the maximum reflected voltage amplitudes in the simulation are approximate 
to the values in the theoretical, and the time value in simulation is shorter due to the decrease of the velocity 
in propagation [27]. The comparison of the times to receive the maximum reflected voltage amplitudes 
between simulation and theoretical calculation with three structures is in Table 4. 
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Amplitude(V) 
o 


Time(s) x 10% 
Figure 6. The amplitude of input and in the FEM simulation 


Table 4. Comparison the simulation results with the theories for time to reach the largest amplitude of 
reflected waves in three structures 


Parameters SAWI1 SAW2 SAW3 
az (um) 870 939.6 1009.2 
a3 (um) 452.4 522 591.6 

Ti (Us) 0.487 0.526 0.596 
Tıı(us) 0.5157 0.5533 0.5909 
T2 (us) 0.806 0.8928 1.04 

T22 (us) 0.798 0.8732 0.9942 
T3 (us) 1.132 1.268 1.475 
T33 (us) 1.0542 1.1673 1.2813 


With the results shown in Table 4, the difference in time to reach the maximum amplitude of response 
waves between the simulation and theory at the third reflector is appropriate because the waves propagating 
from the IDT to the third reflector were weaker and the velocity of the waves decreaseed. The difference in 
distances of these three structures simulated describes the change of material structure in the measurement of 
physical quantities with the distance between the IDT and the first reflector, and between these reflectors in the 
first structure are the smallest, respectively in the third structure are the largest. Table 5 shows results obtained 
from the FEM simulation, and these results are used to input model of the reader in MATLAB, where Gain;, Ti, 
and l; are the amplitude loss, the time delay, and the length of the reflected waves, respectively of the structures 
of SAW1, SAW2, SAW3. Figure 7 shows the signals received from three reflectors of the first structure. The 
output of SAW is a copy sequence of attenuated and delayed from the input. 


Table 5. The synthesis of results from the FEM simulation 


Parameter SAWI1 SAW2 SAW3 
Gain, (V) 0.05604 0.05569 0.05708 
Gain: (V) 0.02687 0.009257 0.03264 
Gain; (V) 0.01811 0.03811 0.05088 
T (us) 0.487 0.526 0.596 
Tə (us) 0.806 0.8928 1.04 
T3 (us) 1.132 1.268 1.475 
lı (us) 0.236 0.25 0.2112 
1, (us) 0.02236 0.13 0.1369 
l; (us) 0.126 0.065 0.175 


The product of (9) and (10) is: 
Sim(t) = A * cos(w * t) * iL, A* A; * cos (w* (t—7;) + 8) *[](t —7,) (11) 
in (10) corresponding to the input signal of analog filter 1, described in Figure 8(a). 
Shifting the phase of (9) and (10) by 90 degrees, then multiplying them together (corresponding to 


the input signal of analog filter2), we get: 


SQm(t) = A * cos( w * t — 7) * 0}, A * Aj *cos(w * (t — Ai) + 0) [[(t- t) (12) 
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Figure 8(b) shows the graph of Som (t). From formula: 


cos(x) *cos(y) = ; [cos(x +y) + cos(x — y)] 


In (11) becomes: 


Sim(t) = A? * A; * Y2, T(t — Ti) * 5x [ cos (w *t +w * (t— T;) +8) + cos (wt; + 8)](13) 
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Figure 7. Response signals from 3 reflectors of the SAW1 
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and (12) becomes: 

Som(t) = A? * Aj * Dey [Ct -— 7) * ~*[ cos (@*t— Sto (t= T)+0+ cos (wT; +0)] (14) 
after low pass filtering Sym (t) and Som (t) we get the signals in Figures 9(a) and 9(b), respectively with (15), 
and (16) where the phase shift can be calculated by sampling the signal levels of in-phase and quadrature 
channels. 

Sm (t) = Lier Ci * [ICt— Ti) * cos (wt; + 0) (15) 


Som(t) = Dizi Ci * [I(t — T:) * sin (wr; + 0) (16) 


Where C; = > * A? * Aj. 
From (15), (16), Sim(t) and Sgm(t) consist of separate segment independent sine and cosine signals, 
we consider any component corresponding with the coefficient i. 


Simi = Ci * COS( w T; + 8) (17) 


Somi = Ci * sin(w T; + 0) (18) 
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Figure 9. Low-pass filter in-phase quadrature of SAW1 for (a) in-phase and (b) quadrature 
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From (17) and (18) we get: 


CG = [szi + Sómi (19) 


where C; is the coefficient of reflectors in output. The output signal amplitudes from three reflectors are 
described in Figure 10. 
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Figure 10. The output signal amplitudes of three reflectors 


With the parameters of structure | in Table 4, let the phase delay angle of the reflector 1, 2, 3 be 91 
P2, and p3 respectively. From MATLAB/Simulink we get sin(@;) and cos(;). Because @; is larger than 
2m, arccos(q@;) will give the incorrect value of @;. In this paper, the distance of the first-to-second reflector 
(dz) is equal to the distance of the second- to-third reflector (d32), we get T21 = Tz — T1 and T32 = T3 — Tp, 
where T31 is the time delay between reflector | and reflector 2, T32 is the time delay between reflector 2 and 
reflector 3. Because d32 = dz,, when the temperature or physical parameters change, the difference g32 = 
21 is very small and it is in the range of 2m which calculates arccos(@32 — #21) accurately. We have 
the (20): 


sin[(p3 — P2) — (P2 - 91)] 
= sin(Y3 — P2) * COS (P2 — P1) — COS(P3 — P2) * sin (P2 — P1) 
= [ sin(p3) * (cos(2))* * cos (p1) — cos (p3) * sin (p2) * cos(p2) 
* COS(Y1) + sin(3) * cos(Pz) * sin (p2) * sin (p1) — cos(P3) (20) 
« (sin(~2))* * sin (9,)] 
— [cos(p3) * cos(pz) * sin (p2) * cos(~1) + sin(ps) * (sin(y2))* 
* CoS(P1) — cos(3) * (cos(~z))* * sin(p1) — sin(ps) * sin(p2) 
* COS(P2) * sin (P;)] 


where sin(@,), cos (@,), sin(@z) cos(@z) sin(@z), cos(@3) are obtained from the model on MATLAB, 


we get sin[(~3 — p2) — (p2 — p1 )] = 0.9754 so [(p3 — 2) — (G2 — Y1)] = 4.49. Similarly, we have 
phase values for SAW2 and SAW3. Table 6 shows the phase shift angles of three structures. 


Table 6. The phase shift angles of three structures of SAW 
Sensor [(@3-2) - (p2 -1)] 


SAW1 4.49 
SAW2 5.39 
SAW3 5.78 
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Table 6 shows that as the distance between the IDT and the reflectors increases, the phase shift 
angles also increase. The effect of physical factors such as temperature, stress, and pressure on the sensor 
cause changes in the propagation length and velocity of SAW, which change the time delay and the phase 
angle in the reflection [24]. 


4. CONCLUSION 

In this paper, we combine the parameters from the SAW model of FEM simulation with blocks of 
the SAW model in MATLABwith three structures corresponding to three different positions of the reflectors, 
which shows the influence of the measurands on the sensor. To eliminate the ambiguity when the phase shift 
exceeds 360 degrees, we used the time delay 132, between the 3"! and 2" reflectors, and the delay time, 121, 
between the 2nd and Ist reflectors. We propose an algorithm to calculate the phase shift between the 
reflected pulses, with the increasing distance of the position between the IDT and the first reflector as 870 
micrometers, 939.6 micrometers, and 1009.2 micrometers, corresponding to phase angles 4.49 degrees, 5.39 
degrees, and 5.78 degrees also increased gradually. The combination of simulation on two softwares is of 
great significance in manufacturing and testing the reader of the passive wireless SAW sensors. 
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